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We study the spatial distributions of two randomly interacting species, in the presence of 
an external multiplicative colored noise. The dynamics of the ecosystem is described by 
a coupled map lattice model. We find a nonmonotonic behavior in the formation of large 
scale spatial correlations as a function of the multiplicative colored noise intensity. This 
behavior is shifted towards higher values of the noise intensity for increasing correlation 
time of the noise. 
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1. Introduction 

Recently there is an increasing research activity on the role of noise in spatially 
extended systems PEE]- One of the challenging problems to be understood is the 
formation of large scale spatial correlations. In this work we present a theoretical 
investigation on the role played by the colored noise in the transient dynamics of 
the spatial distributions of two competing species. The colored noise in spatially 
extended systems has been investigated in the past (and references cited there) , 
0]. Aim of the work is to present a stochastic model with a realistic noise and 
to suggest a possibility of application to biological systems, where the presence of 
fluctuations, such as random variability of temperature, can modify strongly the 
dynamics of an ecosystem 000 . The theoretical results in fact may be compared 
with the experimental data obtained for the biomass of species measured in the Sicily 
Channel [S] . Our investigation could contribute (i) to select environmental forcings 
which affect the dynamics of biological systems 6,9,10 , and (ii) to implement 
physical models for the interpretation of spatial distributions of marine species. 
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2. The Model 

To study the colored noise effect on the formation of large scale spatial correlations 
of two species densities, we consider a coupled map lattice model |11II12| 



7^(1 - - P n yf tj ) + x^j) x ~ (1) 



8 

where x™j and yfj are the densities of two species in the site at the time step 
n, 7 is proportional to the growth rate, D is the diffusion coefficient, Y2s * s * ne sum 
over the four nearest neighbors, and (3 n is the interaction parameter. Here {C,fj)x 
and (Ci'j)y are Ornstein-Uhlenbeck colored noises at each site and correspond 

to the value of the stochastic continuous process Q (t) given by the equation 

^ = --C/(t) + -&(*), (! = *,V) (3) 
at t c t c 

taken at the time step n. Here £j(t) (i = x,y) are Gaussian white noises within 
the Ito scheme with zero mean and correlation function (£,i(t)£,k(t')) = <?5(t — t')8ik- 
The correlation function of the processes of Eq. © is 

<0(*)a(*')> = ^e-l*-*'!/^ (4) 
2r c 

and gives cr6(t — t')6ik in the limit r c — > 0. In absence of noise, for (3 < 1 both 
species survives in a coexistence regime, while for (3 > 1 one of the two species 
vanishes after a certain time and an exclusion regime takes places. To describe the 
noisy interaction between species and the environment we consider an Ito stochastic 
differential equation for the parameter [3 with a periodical driving force jl2lllHIII4] 

d/3(t) dU{(3) . , _ , , 

dt = d(3 + acos ^ + ( 5 ) 

Here U (/?) describes the two stable states of the ecological system (coexistence and 
exclusion) 

U(j3) = h((3 ~ (1 + p)fl^ - 2h(/3 - (1 + p)) 2 /^ 2 . (6) 

The values of the parameters h, p and ?y depend on intrinsic biological characteristics 
of the species, as stability of their life cycles respect to external perturbations and 
ability of one species to compensate the increase or decrease of the other species. 
The periodical driving mimics the climatic temperature oscillations, and £p(t) is a 
white Gaussian noise with < £p(t) > = and < £,3(i)£/3(i') >= vp6(t — The 
interaction parameter (3 n of Eqs. QJ, corresponds to the value of the stochastic 
process (3{t) of Eq. (JSJ) taken at the time step n. In the absence of the additive noise, 
(3(t) oscillates periodically and its values always remain below (3 = 1, giving rise to a 
coexistence regime. We fix the additive noise intensity at the value ap = 2.65 • 10~ 3 
corresponding to a competition regime with (3 periodically switching from coexis- 
tence to exclusion regions. This produces the appearance of stochastic resonance 
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(SR) effect in the dynamics of interaction parameter and noise-induced anticorre- 
lated periodic oscillations in the time evolution of the two species densities |12II14| . 
We analyze the spatial distribution of the ecosystem in this SR dynamical regime, 
by varying both the intensities of multiplicative colored noise in Eqs. (1) and (2). 

2.1. Spatial distributions 

By using a square spatial lattice with N = 100 x 100 sites, and species densities 
vanishing on the grid borders, we obtain the spatio-temporal distributions of the 
two species for uniform initial condition and for three values of the correlation time 
r c , namely r c = 10~ 2 , 10, 10 2 . We find that for the lowest value of the correlation 
time (t c = 10~ 2 ), the spatial distribution is almost the same as that found for 
multiplicative white noise In Fig.^the spatial distributions of the two species 




Fig 1. Spatial distributions with r c = 10 2 at different times for (a) a = 10 12 , (b) a = 10 4 , (c) 
a = 10 5 . The values of the parameter set are: erg = 2.65 • 10 -3 , 7 = 2, D = 0.05, a = 1.5 ■ 10 _1 , 
lu /(2tt) = 10~ 2 , /3(0) = 0.94, N = 100 X 100, and x° . = y9 = 0.5 at all sites (i,j). 

are obtained with correlation time r c = 10 2 and for a — 10~ 12 , 10~ 4 , 10 5 . We see 
that for very low noise intensity (Fig.QJa)) the spatial distributions of the species 
are uniform as at t = 0. For higher noise intensities (see Fig. [IJb)) an anticorrela- 
tion between the two species is observed: the two species tend to occupy different 
positions. Further increase of the noise causes the anticorrelation to disappear and 
the two species densities become uncorrelated (see Fig. dc)). We observe that an- 
ticorrelated behavior is stronger for increasing r c . In particular by comparing the 
distributions obtained with different values of correlation time t c we note that for 
t c = 10 2 the two species densities are strongly anticorrelated. This is expressed by 
the presence, in the same spatial region, of a high density (dark zones) of species 1 
(species 2) corresponding to a low density (light zones) of species 2 (species 1) (see 
Figs. [Hb)). 
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2.2. Correlation between species 

To analyze the role of the colored noise on the correlation between the two species, 
we calculate for different values of the correlation time r c , at the time step n, the 
correlation coefficient c n defined on the lattice as 



(7) 



with 



N 



(8) 



where x n , s™, y n , s™ are the mean value and the root mean square respectively of 
species 1 and species 2, obtained over the whole spatial grid at the time step n, 
cov™ y is the corresponding covariance and N = 100 x 100 is the number of sites of 
the considered lattice. The behavior of the correlation coefficient c n as a function 
of time has been obtained for r c = 10~ 2 ,10,10 2 , and for different levels of the 
multiplicative colored noise. For a = the species are strongly correlated and c™ 
is constant. In Fig. [21 we report the behavior of c™ for r c = 10~ 2 . For low levels of 
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Fig 2. (a) Correlation coefficient c n as a function of time, for r c = 10 2 and for different values of 



the multiplicative noise intensity <r, namely a 



10" 



,10" 



, 10" 



, 10" 4 , 10+ 1 , 10+ 3 ; (b) time 



average correlation coefficient < c n >t as a function of the multiplicative noise intensity. 

the multiplicative colored noise (a = 10~ 14 , a = 10- 11 ), the values of c™ oscillate 
between 0.5 and 0.8, showing a stronger correlation between species for lower values 
of a. For intermediate values of noise intensity, namely for a = 10~ 10 ,10~ 4 the 
spatial correlation coefficient c n oscillates between —0.6 and —0.2 at the frequency 
of the periodical forcing, showing strong anticorrelation between the two species. By 
increasing the noise intensity, the anticorrelation is reduced and finally it disappears 
(cr = 10 +1 , fj = 10 +3 ), that is the species are totally uncorrelated. A nonmonotonic 
behavior of the correlation coefficient as a function of the multiplicative intensity 
appears clearly. This behavior is more evident by considering the time average 
correlation coefficient < c" > t as a function of the multiplicative noise intensity (see 
Fig- Gib)). The minimum in Fig.E^b) corresponds to the anticorrelated oscillations 
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shown in Fig. HJa). Moreover this behavior is connected with the anticorrelated 
oscillations present in the time evolution of two competing species in each point of 
our spatial grid |14j . By comparing the results shown in Fig. 2 with the white noise 
case we observe that the effect of a small correlation time r c is to produce an 
enlargement of the minimum of < c™ > t towards higher values of the multiplicative 
noise intensity. The ecosystem remains in his anticorrelated regime for a longer time 
due to the memory of the colored noise. The two noise sources have different roles 
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Fig 3. (a) Correlation coefficient c n as a function of time, for r c = 10 2 and for different values 
of the multiplicative noise intensity a, namely a = 10 -11 , 10 -8 , 10 -5 , 10 +1 , 10+ 3 , 10+ 5 ; (b) time 
average correlation coefficient < c™ >t as a function of the multiplicative noise intensity. 



on the dynamics of our ecosystem. Specifically the additive noise affects directly the 
dynamical regime, while the multiplicative noise breaks the symmetric dynamical 
behavior of the two species by producing a coherent response of the system [12111 1| , 
which is responsible for the appearance of anticorrelation in the spatial distributions 
of the two species. In Fig.|3the behavior of c™ and < c™ > t is reported for t c = 10 2 . 
As in the previous case we observe a nonmonotonic behavior of c" as a function of 
the multiplicative noise intensity. For low noise intensities (a = 10" 11 , a = 10 -10 ), 
the species are strongly correlated and c™ is almost constant. For a — 10~ 8 strong 
oscillations between +1.0 and —1.0 appear. This effect is due to the colored noise, 
in fact the ecosystem remains periodically in the correlated or anticorrelated regime 
for a time interval of the order of the correlation time of the noise. For greater noise 
intensity (a = 10~ 5 , 10 1 ), the anticorrelation regime takes place (c™ w —1.0). By 
increasing the noise intensity (er = 10 +3 ,10 +5 ), c™ goes through a nonmonotonic 
behavior, taking on values around zero (uncorrelated regime). Comparing Fig. [21 
and Fig. we observe that the presence of a bigger correlation time r c determines 
a stronger anticorrelation in the two species distributions. The minimum of the 
time average correlation time < c™ > t is deeper and it is shifted towards higher 
values of the multiplicative noise intensity for greater values of correlation time r c . 
A highly correlated noise therefore produces two effects: (i) the system is pushed 
towards a strongly anticorrelated regime (c™ « —1-0); (ii) this regime takes place for 
higher values of multiplicative noise intensity. The colored noise produces therefore 
a persistence of the dynamical regime experienced by the ecosystem. 
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3. Conclusions 

In this work we analyze the effects of the colored noise in the spatial distributions 
of two competing species. By using a discrete time evolution model, which is the 
discrete version of the Lotka-Volterra equations with diffusive terms, in the pres- 
ence of a multiplicative colored noise, with a random interaction parameter, we 
analyze the spatio-temporal behavior of the two species. We find that the colored 
noise determines an enhancement of the effects observed in the presence of white 
noise ^2] > that is (i) formation of large scale spatial correlations with the same peri- 
odicity of the deterministic force and (ii) anticorrelation in the spatial distributions 
of the two species. Finally we observe that (i) the intensity of the anticorrelation 
increases as a function of r c , (ii) the system remains in the anticorrelated region for 
a bigger range of multiplicative noise values as r c increases. The noise spectrum 
of a real physical system is characterized by a cut-off (correlated noise). Therefore 
data from real ecosystems, whose dynamics is strongly affected by fluctuations and 
by environmental physical variables, can be modeled using sources of colored noise. 
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